Nonequilibrium tricriticality in one dimension 
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We show the existence of a nonequilibrium tricritical point induced by a repulsive interaction 
in one dimensional asymmetric exclusion processes. The tricritical point is associated with the 
particle-hole symmetry breaking introduced by the repulsion. The phase diagram and the crossover 
in the neighbourhood of the tricritical point for the shock formation at one of the boundaries are 
determined. 



The lack of any free energy like entity to describe 
nonequilibrium steady states makes general studies of 
nonequilibrium phase transitions rather difficult. For this 
reason, cvcnthough different types of phase diagrams are 
known with first-order or continuous phase transitions, 
and critical points from various case studies (EH, SflU, 
generalizations of these results are not straightforward. 
Since steady states are the closest analogs of equilibrium 
states (both being stationary in time), it is important 
to know how far the richness of the equilibrium critical 
phenomena with their connections to symmetries Q can 
be found in nonequilibrium systems. 

This paper shows the bifurcation of a line of critical 
points through a tricritical point 0, 0] in the phase dia- 
gram of a one dimensional interacting driven system as 
one interaction parameter is tuned. A tricritical point is 
also a point where a first order transition changes over 
to a continuous one, the pair of critical lines being the 
edges of two first order surfaces. We do observe such 
equilibrium-like gross features in the nonequilibrium sys- 
tem so as to classify it as a tricritical point, but, unlike 
equilibrium systems, all happening in one dimensional^ . 

The occurrence of a tricritical point in any system is 
significant because it implies a confluence of two different 
phenomena. As a consequence, the tricritical point con- 
trols the scaling description in its neighbourhood. Known 
examples from equilibrium systems include the #-point 
of polymer solutions |9|, bunching and phase separation 
of steps on crystal surfaces [lfj, different transitions in 
magnets @, current hunt for a tricritical point in quark- 
gluon plasma in the context of early universe 11( , and 
many othersQ. In equilibrium systems, a tricritical point 
is one in the hierarchy of critical points and corresponds 
to the case of three relevant variables as opposed to two 
for an ordinary critical point, as e.g., temperature and 
magnetic field for a Curie point of a magnet. In the Lan- 
dau theory of phase transitions, the sequence of critical, 
tricritical points ... occur as new symmetric minima in 
the Landau function develop. In analogy with that, we 
find in the nonequilibrium problem in hand, the tricriti- 
cality is associated with a particle-hole symmetry break- 
ing or successive occurrences of peaks in the current in 
the system, that allows maximum current through the 
system at two different densities. 

Our results are based on a class of well-studied one- 
dimensional models for interacting particles moving on 
a track, which arc variants of the asymmetric exclusion 



process (ASEP). There is current interest in these mod- 
els because of the recent observation of localized shock 
phases and associated criticality. A shock in these mod- 
els is a discontinuity in the steady state density along the 
track. Let us consider a onc-dimcnsional lattice with par- 
ticles injected at site i = at a rate a and withdrawn at 
i = N at a rate 1 — 7. The particles hop to the right with 
mutual exclusion so that the occupation number at a site 
is or 1. In addition, there is a next nearest neighbour 
repulsion[l2j so that a configuration 
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with a rate 1 
at a rate 1 — 1 



with (0 < e < 1), where occupied and unoccupied states 
of a site are represented by 1 and respectively. The 
hopping process is particle conserving. With a nonvan- 
ishing current, the system can evolve to a nonequilibrium 
steady state. Nonconservation is introduced by allowing 
evaporation or desorption (1 — > 0) at a rate u>d and de- 
position or adsorption (0 — > 1) at a rate uj a of particles 
at any site on the track. This dynamics, called the Lang- 
muir dynamics, can maintain a density 



UJa , 1 ., , 

PL = ; + ~ if U a ^ UJ d . 



(1) 



This is an equilibrium-like process with no current in 
the system. If the net flux due to adsorption/desorption 
is comparable to the hopping current, there will be a 
competition between the attempt to equilibrate and the 
drive to the nonequilibrium steady state. As a result, one 
gets new features like a shock phase, critical points one 
would not get otherwise. 




FIG. 1: (a) The critical lines in the 7-r plane, where r = e— ej. 
The width vanishes with a characteristic exponent as r — » 0—. 
(b) Three dimensional view of the tricritical point J in the a- 
7-r space. Only the critical lines are shown. These lines are 
nonplanar and (a) shows the projection. 
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For large N, we may use a continuum notation x = i/N 
and the average occupation of a site becomes the density 
p{x). For given microscopic parameters for the hopping 
rules and nonconservation, the steady-state density pro- 
file p(x) depends on the external parameters a and 7. 
The possible noncquilibrium phases are then represented 
in a phase diagram in the a-7-e space. 

The extra interaction e is important in particle-hole 
symmetry breaking. To see this let us consider the con- 
served case (u) a = ujd = 0). Instead of particle injection 
at x = 0, we may consider holes being injected at the 
right (x = 1) at a rate 1 — 7, hopping left with iden- 
tical rates and withdrawn at x = at a rate a. This 
particle-hole symmetry allows a symmetric form for cur- 
rent j(p) = j(l — p). For e = 0, j(p) = (1 — p)p satisfies 
the symmetry with a maximum at p = 1/2. A similar 
single peaked current with j(p) = j(l — p) is expected for 
low interaction strengths (small e). If a system evolves 
to the maximal current state then the steady state or the 
phase also has the particle- hole symmetry (0 <-> 1). In 
addition to this symmetric phase, there are two distinct 
phases related by the <-> 1 symmetry but separated in 
the phase diagram by a first order line. In the other ex- 
treme case £ = 1 which forbids the hops 0101 — > 0011, 
there is a zero current state for p = 1/2 as for p = 
(empty track) and p = 1 (packed track). Therefore, for 
strong enough interaction, p = 1/2 is a local minimum 
of the current while the maximal current carrying states 
are symmetrically located on the two sides of p = 1/2. 
The exactly known stationary current for the conserved 
system does show the change from single to double peak 
for e close to e = £j ~ 0.8 [13, EH- In any one of the 
maximal current state, the particle-hole symmetry is not 
respected [l9|. The overall symmetry is recovered by rec- 
ognizing that the <-> 1 transformation gives the other 
maximal current state. The single maximal current phase 
in the single-peak case now breaks up into several phases, 
two of which are the two maximal current phases. The 
occurrence of the two maximal current phases in differ- 
ent regions of the phase diagram is a reflection of the 
symmetry breaking introduced by the interaction. 

For the ASEP case (i.e., with e = in the above 
model) @, H, 3 , for w a 7^ <^->d, there are three types 
of phases, namely, (i) an injection rate controlled, to 
be called the a-phase, (ii) a withdrawal rate controlled 
phase, to be called the 7-phase, and, most importantly, 
(iii) a phase containing a localized shock, to be called 
a shock phase. The maximal current phase of the con- 
served case becomes unstable under nonconservation, ex- 
cept for the case of accidental symmetry at uj a = uJd- The 
major signature of nonconservation is in the appearance 
of the shock phase with a localized shock in the density 
profile in between the a- and the 7-phases. This phase 
replaces the first-order boundary of the conserved case. 
The maximal current density is still important because 
the shock is centered around the maximal current density. 
In other words, a shock connects configurations related 
by the particle hole symmetry, i.e., density profiles with 



p = 1/2 ± 5. The transition to the shock phase is gener- 
ally first order because the shock height is nonzero at the 
transition, but critical points do occur where the shock 
height vanishes. The vanishing height implies a special 
state that regains the particle-hole symmetry. Such a 
critical point therefore has to be at the peak of the cur- 
rent. There are characteristic universal exponents, for 
both bulk [2| and boundary Q, associated with the criti- 
cal point. 

In the interacting case (e > 0) with conservation (uj a = 
u>d = 0), a very complex phase diagram with seven differ- 
ent phases including two maximal current phases (men- 
tioned earlier) is known but without any shock phase. 
let alone the critical point. The dynamics of the inter- 
acting model has also been studied in certain regions of 
the parameter space in Refs. [H, EH- With nonconser- 
vation, localized shocks appear, and double shocks and 
downward shocks have also been seen[3, EH, EH but the 
phase diagram is not known. In the broken symmetry 
case, a shock centered around one peak will not show the 
<-> 1 symmetry and in fact the state obtained by such 
transformation would occur in a different region of the 
parameter space. Given the fact that the maximal cur- 
rent densities play an important role in shock formation, 
the symmetry-breaking is expected to show new critical 
points compared to the symmetric case. Consequently, 
the critical point for the shock phase for e < £j appears 
as a line in the extended a-7-£ phase diagram and this 
line undergoes a bifurcation at £ = £j where the symme- 
try breaking takes place. The bifurcation point J is the 
tricritical point as shown in Fig. [TJ and it has its own 
distinct scaling. 

The phase transitions can be analyzed by using a 
boundary layer analysis for the density profile in the con- 
tinuum, long time, long length scale limit, the so-called 
hydrodynamic regime [J, Il7j . The shock phase can be 
seen as forming via a thickening and eventual decon- 
finement of a boundary layer ('shockening transition'). 
Moreover, this bulk cum boundary shockening transition 
is in turn associated with a dual boundary transition 
where the boundary layer changes from a depletion region 
to an accumulated region 0, E3 ■ If the shockening and 
the dual transition lines intersect then there is a critical 
point ( "self-dual" ) . 

So far as the single to double peak change in the cur- 
rent is concerned, a Taylor series expansion in p and £ of 
the exactly known current [e| upto fourth order is suf- 
ficient. For simplicity we work with this expansion in a 
general form 

■1 \ 2 ?' + u r. 1 2 l- 4 
3(p) = — -- 2 (p-- 2 ) 2 -u(p--)\ (2) 

with the constant piece chosen to ensure that the current 
vanishes for p = 0, 1 (the empty track and the fully- 
occupied track). Eq. J2J for r = 2,u = recovers the 
known current density for ASEP processes (J = (1 — p)p) 
while the double peak appears for r < 0. Throughout u 
is taken as a positive constant, and r is tunable but kept 
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FIG. 2: (a) Current versus p for three values of r. (b) shows 
the notations used for r < 0. The peak densities are p c and 
1 — /9 C - The inner layer connects pairs of densities (with same 
current) from p ,pi, 1 — pi,l — p - For pi = 1/2, the low 
and high densities are p* and 1 — p* . For currents les than 
j{l/2) (and also for r > 0), p a and 1 — p are the two relevant 
densities. 



small. The current is shown in Fig. [2^ for typical values 
of r. In the double peak case, there are a few special 
densities we need, namely, (i) the densities for maximal 
current, p c , l—p c , and (ii) the densities p* , 1— p* , with the 
same current as the minimal one (j(l/2)). Any constant 
current line with j(p c ) > ] > j(l/2) intersects at four 
densities p a , pi, 1 — pi, 1 — p a . There are only two values 
p , 1 — p if j < j (1/2) ■ These are shown in Fig. [2t>- 

The hydrodynamic approach based on the continuity 
equation with j (p) of the conserved case is known to agree 
remarkably with numerical simulations (l2l [Til ]. In this 
approach the density satisfies a continuity equation sup- 
plemented by an additional contribution from the non- 
conserving process. The density then satisfies 



dp 

at 



dJ_ 

dx 



(3) 



where J is the current consisting of the bulk current j(p) 
and a diffusive current, and 



So(p) = — fi(/0L - p), ft = (w + Ud)N, 



(4) 



accounts for the Langmuir dynamics. f2 is related to the 
net flux of particles and pl is the density defined in Eq. 
[TJ Nonconservation would matter only when the net flux 
of particles adsorbed or desorbed is comparable to the 
current in the system. For this reason CI is kept constant, 
finite (scaling limit) for N — ► oo. By adding the diffusive 
(Fick's law) current, J in Eq. [3] can be written as 



J 



(5) 



with the diffusion constant e ~ 0(1/ N) is a small param- 
eter (not to be confused with e), and is a remnant of the 
underlying lattice. This form, Eq. O also follows from a 
continuum limit of the discrete model. Despite its small- 
ness e plays an important role in the phase transitions, 
especially in shock formation. The steady state density 
profile is now given by 



' dx 



So(p) = 0, 



where Si(p) 



dj(p) 
dp ' 



(6a) 
(6b) 



and the boundary conditions p(0) = a, and p(l) = 7. 
We already noted that p = 1/2 is a special density due to 
the particle-hole symmetry. To avoid extra complications 
arising from the accidental matching of densities, we take 
uj a LOd so that pl 7^ 1/2. Such cases will be considered 
elsewhere. 

Eq. [5a] entails two length scales, (i) x for the bulk and 
(ii) x = (x — xq)/e which is significant in a thin region as 
e — > around an appropriately chosen x$. Eq. [6al does 
not depend on the choice of xq . The separation of the two 
scales is used to develop a uniform approximation of the 
solution order by order in e. Here we restrict ourselves 
to the lowest order approximant. The bulk solution in 
terms of x comes from the first order equation obtained 
with e = in Eq. [Sa| This solution, to be called the 
outer solution, 



p(x) = Pout (a;) 



is given implicitly by 



g r (p) = 2ft x + C, 



(7) 



(8a) 



where C is a constant to be fixed by only one of the two 
boundary conditions, and 



9r{p) 



(uY*+r)Y+^Y h Y- 2 + ^Y 3 
+ (uY* +r)Y L \n\Y h -Y 



(8b) 



with Y = 2p— l,and Yl = %Pl — 1- We shall suppress the 
r-dependence in g r , unless required. For the left bound- 
ary condition, the outer solution is 



with a density 



determined by 



g(p) = 2nx + g(a), 



Po = Pout(l) at x = 1, 



g( Po ) = 2Q + g(a). 



(9) 



(10) 



(11) 



Here, p D depends on a, but in general, p D 7^ 7. 

To satisfy the other boundary condition p = 7 at x — 
1, we use the second scale x = (x — l)/e around the 
boundary point x = 1. With this variable, the density 
(to be called the inner solution) satisfies 



d 2 p in 
dx 2 



Si(p-n 



) d£in 

dx 



0, 



(12) 



obtained from Eq. (f6"a| by first changing the variable to 
x and then taking e — ► 0. The absence of the Sq term can 
be understood from the observation that for a constant f2 
this layer is too thin for the nonconservation to matter, 
at least to leading order in e. For a smooth density profile 
(the solution of the second order equation with any e > 
no matter how small), we need to match the outer and 
the inner solutions by requiring that 



lim Pout (a) 

x— ►! 



lim p in (£). 



(13) 
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FIG. 3: Schematic diagram of the phase boundary for shock 
at x = 1. Thin, thick and dotted lines represent continuous 
transitions, first order transitions and dual lines respectively. 
"A"["D"]: accumulated [depleted] bounary layer. C, Ci, C2 
are the critical points and J the tricritical point, "a" and 
"s" represent a-phase and shock phase with 1,2 denoting 1- 
shock, 2-shock phases. The phase boundary separating the 
a and the Is phases is given by Eq. (|20|l . Downward shocks 
appear in the hatched region (with 7 < 1/2) in (c). The dual 
lines demarcating "A"/"D" are not shown in (c). 

Here x — > — 00 gives the outer limit of the inner solution. 
Incorporating this matching condition, Eq. [T2] can be 
written as 

^T=j(Pin)-j(Po), (14) 

in the thin boundary layer, with p- m (x = 0) = 7. The 
complete matched solution is obtained by joining the in- 
ner and outer approximations and subtracting their com- 
mon value. Therefore, the density profile is given by 

p(x) = p ut (a:) - Po + Pin(x) + 0(e). (15) 

Eq. [15] identifies the scale dependent separation of the 
bulk and the boundary contributions and it provides a 
uniform approximation of the density in the whole do- 
main including the boundaries. 

Let us first consider the r > case. The current has 
only one peak at p = 1/2. The transport across the track 
is analogous to the well-understood ASEP case. For a 
low injection rate a, the bulk density profile does not 
depend significantly on the withdrawal rate so long 7, 
the accumulated density at x = 1, is small. One gets the 
a-phase. For low 7, there is a depletion layer at x = 1, 
which changes for 7 > 7( j = p (a) to an accumulated 
layer. This is the dual boundary transition mentioned 
earlier. By symmetry (Fig. [5]) and from Eq. (|14[) . p- m (x) 
would saturate at p = \ — p {a) as x — > 00. Consequently, 
for 7 > 7,s = 1 — Po{ol), a boundary layer given by Eq. 
(JU) is not sufficient to satisfy the boundary condition at 
x = 1. In this situation, the density profile is given by 
the two outer solutions, 

g(p Igt )=2Q, ( x -l)+0( 7 ) 

on the right side and 

g(put) = 2£! x + g(a) 

on the left. These two solutions are joined smoothly in 
a thin region by the inner solution of Eq. (|12[) with 



x = (x — x s )/e. On the bulk scale this looks like a discon- 
tinuity and is therefore a shock at x = x s . The matching 
conditions and the location of the shock (x s ) are given 
by 

pm(x — ► x a —) = lira pin(x), (16a) 

x — > — 00 

and p rg t(a; — > x s +) = lira pi a (x), (16b) 

x — >oo 

where Eq. I16al is the matching of the left outer solution 
to the left end of the inner solution while Eq. Il6bl is the 
matching for the right side. This deconfinement of the 
boundary layer by shifting into the bulk as a shock is the 
shockening transition. For a given a, the shock and the 
dual transitions are given by 

7 = 7 8 = 1 - p (a) (shockening) (17) 
7 = Id = Po{a) (dual) (18) 

The critical point (a c , 7c ), being the point of intersection 
of the two lines, has j c — 1/2, and a c satisfying Eq. (fTTj) 
with p = 7 C . 

The shape of the transition curve near (cv c ,7c) for any 
r > follows from Eq. (fTTj) , by noting that dg(p)/dp is 
zero at p = 7 C but not at p = a c . This shows that the 
shape exponent x is given byQ 

7 -7c H a-a c | x , with x = 1/2, (19) 

close to the critical point. The phase diagram is shown 
in Fig. [3^,. The critical point (a c (r), 7c ) changes with r 
but the exponents remain the same as in the ASEP case. 
We note in passing that for 7 < 7c , one crosses the dual 
phase boundary so that there is a depletion layer near 
x = 1. This layer is non-shockening in the sense that 
it never reaches saturation 0, [T3- It acts as a shield for 
the boundary density, rendering the "effective" boundary 
condition at x = 1 as 7 = 7c . Consequently, thanks to 
this depletion layer, there is a continuation of the critical 
point in the phase diagram. 

The situation is completely different for r < 0. We 
refer to Fig. Let us start with small injection rates. 
There is a dual line at 7 = p (ct) so long p (a) < p* 
(Fig. [2]). There will be a symmetric shock centered 
around p — 1/2. The transition to the shock phase (or 
the shockening transition) is at 7 = 7 S = 1 —p (a). How- 
ever, with change of a when p a exceeds p* , the boundary 
layer bridges p to p\ < 1/2. The dual line continues to 
be 7 — "fd — p (a), but the shock line is given by 

7 = 7s 1 = P\{a). (20) 

The shock is therefore smaller in height, lying below 
p = 1/2. A zero height shock at the boundary is now 
possible if the following two conditions are met: (i) the 
injection rate a is such that the outer solution gives the 
peak density at the boundary, i.e., p D (a) = p c , and (ii) 
the withdrawal rate is such that the accumulated density 
7 is also p c . Under these conditions, the density will then 
have an infinite slope at the x = 1 boundary. Therefore, 
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7 — Pc — Po (a) gives a critical point. This locates an 
off-center critical point Ci in Fig ^p, with 

7d = Pc and a = a cl , (21) 

where a c i satisfies Eq. [IT] with p Q = j c i. 

The density profile close to this critical point C\ (for 
r < 0) is shown in Fig. |4ji. These are obtained by a 
numerical solution of Eq. (j6"a|) in the limit of small e. 
For g > 7 > 7 c i, the shock height starts increasing with 
7 until the shock connects p = p* and p = 1/2. This 
is the largest shock one can produce from the left peak 
of the current. For a range of 7 there is the possibility 
of a downward shock. One sees an overcrowded region 
sandwiched between the a-phase and the 7-phase, even 
though both a and 7 are less than 1/2, but the shock 
is centered around 1/2. The downward shock shifts to- 
wards x = 1 as 7 is changed and it disappears as it hits 
the boundary. At this point the downward shock gets 
converted into a depleted boundary layer. With increase 
in 7, the depleted layer undergoes a boundary transition 
to an accumulated layer leading to a second shock from 
the second peak. There is no obvious symmetry relation 
between the two upward shocks except for their centers. 
The sequence is shown in Fig. [4^,. For 7 > 1 — p*, the 
two shocks merge into one bigger one centered around 
p=l/2. 

Next, consider a > a c \ as shown in Fig. With the 
largest shock from the low density peak of the current, 
it is possible for the density to reach the second peak of 
the current, p = 1 — p c , at x = 1. This gives the second 
critical point because the shock height starts from zero 
as 7 exceeds 1 — p c . The second critical point C2 at 
(a C 2,7 C 2) is determined by 

7c2 = 1-Pc, (22a) 
5 ( 7c2 ) = 2ft (l-x.)+g(0), (22b) 
g(0) = 2flx s +g(a C 2), (22c) 




0.9 0.95 X i 0.9 0.95 X 



FIG. 4: Plots of p(x) vs. x (r = -0.2, u = 2.2, fi = 0.1) 
from numerical solution. The values of 7 are indicated on 
the plots while a is kept fixed, close to a c \ in (a) and a C 2 in 
(b). In (a), as 7 increases, one sees a sequence: (i) the shock 
developing from zero height (7 = 0.35) to a maximum (7 = 
0.39), (ii) an additional downward shock (7 = 0.41) which 
goes towards the boundary (7 = 0.56), (iii) the emergence of 
another upward shock of finite height (7 = 0.70) , and (iv) 
finally a bigger shock symmetric around p — 1/2 (7 = 0.71). 
(b) shows the sequence for the second critical point. There is 
no downward shock here, while the first shock is of maximum 
possible height. To be noted is the depletion layer for 7 = 
0.62. 



where x s is the position of the first shock, and g is from 
Eq. © with r < 0. Fig. |3fc shows the nature of the 
phase boundary for a given value of r < 0. These two 
critical points Ci, C2 merge as r — > 0. Given the form of 
g(p) in Eq. f8b|) . we find 

|7ci-7jM7c2-7jI~ M 1/2 , (r-0-), (23) 

where 7j is the critical value at r — 0. The power law 
dependence is expected to be universal. The bifurcation 
at point J is shown in Fig. QJi in the projected plane 
of 7-r. The full three dimensional lines [20| are shown in 
Fig. QJ>. An analysis as done for the r > case shows 
that the shape of the boundaries at these new critical 
points would be similar to that for the r > case, i.e the 
shape exponent will be x = 1/2 as in Eq. (|19|) . implying 
universality. 

Let us now consider the special point J. The physical 
picture developed above remains valid here. The tricrit- 
ical point J, being at the intersection of the shock line 
and the dual line, is at 

5r=o(7j) - 20 + 5r=0 (aj), 7J = 1/2. (24) 

Since Si(p) ~ (p — 1/2) 3 , the shape exponent is 

XJ = 1/4, (25) 

where %j is defined as 

7 — 7j H ot — ckj | XJ . (26) 

The phase boundary is shown in Fig. [3J3. The first order 
transition lines as shown in Fig. [3] become surfaces in 
the three dimensional phase diagram. These first-order 
surfaces end on critical lines and these critical lines are 
shown in Fig. [Tb. The locus of the intersection of the first 
order lines for r < (triangle in Fig [3J;) in the extended 
phase diagram is a line ending at J. This is the line on 
which both the shocks connecting p = p* to p = \/2 and 
p = 1/2 to p = 1 — p* are just at x = 1. This completes 
the identification of J as the tricritical point where a first 
order line is converted into a continuous one. 

The crossover from the tricritical to the critical be- 
haviour is obtained by expanding Eq. [TT] around a = 
O-JiPo = 7,7 an d r — ► 0+. The shape of the transition 
curves and also the height of the shocks near the tricrit- 
ical point can be described by a scaling form 

a-aj ~ |7-7J| 1/XJ T (r 1 7 - 7J I" ' 0J ) , (27) 

where ^"(a;) is a scaling function with the crossover expo- 
nent tpj = 2. For x — > 0, T(x) — > constant while for large 
x, !F(x) ~ x. For r ^ 0, with 7 close to the critical value, 
the scaling variable becomes large and one recovers the 
critical value x _1 = Xj 1 ~ V'J = 2. This also shows the 
importance of the tricritical point because it controls the 
scaling behaviour in its neighbourhood. 

Within the window of the single shock phase for r < 0, 
there is a region of downward shock bounded by critical 
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and first order lines. The peculiarity of the downward 
shock is that it does not traverse the whole lattice but 
goes critical within the track. The thickness of the region 
in the phase diagram is determined by the shallowness of 
the minimum at p = 1/2. Details of these lines and the 
region will be reported elsewhere. 

In summary, we have shown that the repulsive inter- 
action induced particle-hole symmetry breaking in the 
form of a double peaked current in one dimensional asym- 
metric exclusion process leads to a nonequilibrium tri- 
critical point. This tricritical point controls the scaling 
behaviour in its neighbourhood. The tricritical expo- 
nents are different from the critical ones and an addi- 



tional crossover exponent is required for the crossover to 
the critical behaviour. The phase diagram also shows a 
region of downward shock. It remains to be seen if this 
is a necessity for the nonequilibrium tricriticality. 
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